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Abstract —A recent approach in modeling and analysis of the 
supply and demand in heterogeneous wireless cellular networks 
has been the use of two independent Poisson point processes 
(PPPs) for the locations of base stations (BSs) and user equip¬ 
ments (UEs). This popular approach has two major shortcomings. 
First, although the PPP model may be a fitting one for the BS 
locations, it is less adequate for the UE locations mainly due to 
the fact that the model is not adjustable (tunable) to represent the 
severity of the heterogeneity (non-uniformity) in the UE locations. 
Besides, the independence assumption between the two PPPs does 
not capture the often-observed correlation between the UE and 
BS locations. 

This paper presents a novel heterogeneous spatial traffic 
modeling which allows statistical adjustment. Simple and non- 
parameterized, yet sufficiently accurate, measures for capturing 
the traffic characteristics in space are introduced. Only two 
statistical parameters related to the UE distribution, namely, the 
coefficient of variation (the normalized second-moment), of an 
appropriately defined inter-UE distance measure, and correlation 
coefficient (the normalized cross-moment) between UE and BS 
locations, are adjusted to control the degree of heterogeneity 
and the bias towards the BS locations, respectively. This model 
is used in heterogeneous wireless cellular networks (HetNets) 
to demonstrate the impact of heterogeneous and BS-correlated 
traffic on the network performance. This network is called 
HetHetNet since it has two types of heterogeneity: heterogeneity 
in the infrastructure (supply), and heterogeneity in the spatial 
traffic distribution (demand). 

Index Terms —Spatial Traffic Modeling, Heterogeneous Wire¬ 
less Cellular Networks, Stochastic Geometry, Point Processes. 


I. Introduction 

W ITH the advent of the increasingly diversified us¬ 
age scenarios and applications in the envisioned 5G 
wireless networks, the traffic (demand) in time and space 
is getting increasingly heterogeneous (non-uniform). Coping 
with the performance-related challenges in such networks 
often necessitates the availability of realistic, yet relatively 
simple and manageable, traffic models. 
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The statistics of the signal-to-interference-plus-noise ratio 
(SINR) are the key to the performance analysis of hetero¬ 
geneous wireless cellular networks. The signal strengths and 
interference levels depend strongly on the network geometry, 
i.e., the relative positions of the transmitters and the receivers. 
Accordingly, in heterogeneous wireless cellular networks, spa¬ 
tial statistics of the traffic demand (UE distribution) as well as 
those of the service points (BS distribution) have direct impact 
on the network performance. 

In the traditional single-tier homogeneous wireless cellular 
networks, usually, the macro-BSs have been assumed to be 
located in deterministic hexagonal grids (lj. With the advent 
of multi-tier heterogeneous cellular networks (HetNets two 
major transitions have occurred in network modeling. First, 
small-cell BSs with varying communication characteristics 
(such as pico-BSs and femto-BSs) have been envisioned to 
be deployed in the network. It has also become apparent that 
the regular grid topology assumption for BS locations does not 
hold anymore mainly due to the fact that small-cell BSs will 
often be deployed in UE hot-spots which are rather randomly 
distributed. Therefore, in the contemporary HetNets literature, 
BS locations have been modeled by Poisson point processes 
(PPPs). 

In addition to the PPP model for the BS locations, it is 
becoming increasingly common to model the UE locations as 
another independent homogeneous PPP @ 

Although the recently adopted PPP model for the UE 
locations is more realistic in comparison to the most common 
UE location model used in the earlier literature, in which the 
location of a fixed number of UEs in a cell are determined 
through two dimensional (2D) uniform random processes, the 
PPP model does not adequately represent the scenarios in 
which UEs are heterogeneously distributed (e.g., clustered). 
Moreover, the independence assumption between the two PPPs 
(for UE and BS locations) does not capture the correlation 
between the UE and BS locations observed in reality. 

The real UE distributions are seldom pure PPPs. Due to the 
deployment of small-cell BSs in UE hot-spots, the correlation 
between UEs and BSs is an apparent phenomenon in HetNets. 
Network users are usually concentrated at social attractions 
such as residential and office buildings, shopping malls, and 
bus stations. Studying the UE distributions of more extreme 

this paper we use the terms HetNet and Heterogeneous Cellular Network 
(HCN) interchangeably. 
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characteristics and their impact on network performance is 
thus an important investigation. The requirement is a continu¬ 
ously adjustable and tunable traffic model which can represent 
the broad possibilities from completely homogeneous cases 
(e.g., deterministic lattice) to extremely heterogeneous types 
(e.g., highly clustered scenarios), and from BS-independent 
UE locations to highly BS-correlated types. 

Only when this realistic spatial traffic model is used in 
HetNet scenarios, the true impact of traffic characteristics 
on the network performance can be adequately captured. 
Obviously, the spatial UE distribution has a major impact on 
the network’s key performance indicators, such as UE rates 
and outage probabilities which depend directly on UEs’ SINR 
statistics. The spatial UE distribution has also an impact on the 
network energy efficiency. Especially in HetNets with dense 
small-cell deployments, the cell switch-off approach 0-0 
is an effective scheme in energy saving; the extent of cell 
switch-off depends mainly on the distribution of UEs. 

A. Contributions of This Paper 

As a first step to fill the explained void in the literature, the 
contributions of this paper are summarized as follows: 

• A spatial traffic modeling approach with adjustable statis¬ 
tical properties, capturing the severity of the heterogene¬ 
ity and the extent of the correlation with BSs, is intro¬ 
duced. Although the approach is presented in the context 
of heterogeneous wireless cellular networks, it is very 
flexible regarding the underlying network technology and 
it is general enough to be applied to other contexts as 
well, including Wi-Fi, ad-hoc, and sensor networks. 

• Mathematical tools in stochastic geometry, including the 
Voronoi and Delaunay tessellations, are used to illustrate 
the similarities between traffic modeling in the time and 
space domains. This similarity leads to the introduction of 
new geometric inter-point distance measures for capturing 
the properties of spatial point patterns. These measures 
are chosen in such a way that they resemble the well- 
known inter-arrival time (iat) in the time domain. 

• Only two parameters are introduced for a fairly accurate 
description of a heterogeneous and BS-correlated spatial 
traffic scenario: 

1) The coefficient of variation (CoV) values of the 
appropriately defined inter-point distance measures 
are used for specifying the deviations from homo¬ 
geneity. As stated earlier, the discussed measures in 
the 2D space domain can be interpreted as the equiv¬ 
alents of the iat measure in the one dimensional 
(ID) time domain. It is worth mentioning that CoV 
is a commonly used statistic in traffic and queuing 
theories. 

2) The correlation coefficient between the spatial UE 
distribution and the spatial BS distribution is used 
for specifying the bias of UEs towards BSs. 

• The developed methodology is demonstrated in a hetero¬ 
geneous wireless cellular network (HetNet) to illustrate 
the effects of the realistic traffic modeling on the perfor¬ 
mance. 


B. Related Works 

Traffic demand modeling in the time domain has been 
investigated well in the literature 0-0 Traditionally in 
voice-only networks, homogeneous Poissonian models were 
accurate enough to model traffic in time. After the emergence 
of different applications, such as video and data with variable 
rates, the Poisson model failed to capture the traffic statistics 
( 7 ); as a result, various heterogeneous (super-Poisson) traffic 
models based on the hidden Markov model (HMM), Markov 
modulated Poisson process (MMPP) | 8| and other stochastic 
methods have been proposed in the literature and used for 
performance analysis. 

In the space domain, on the other hand, while there are many 
papers concentrating on the modeling of base station locations 
using stochastic geometry © 02 ) there are only few works 
in the literature which take into account the heterogeneous 
spatial distribution of traffic demand in wireless cellular net¬ 
works Q5H22). To the best of the authors’ knowledge, none 
of the existing works provides a statistically adjustable model 
representing a variety of possible scenarios for UE distribution. 

In fl6| , Bettstetter et al., presented an algorithm to cre¬ 
ate a random inhomogeneous node distribution based on a 
neighborhood-dependent thinning approach in a homogeneous 
PPP. The model, however, can not be used for generating BS- 
correlated UE patterns, as this is beyond the scope of that 
model. 

In (17), Qvarfordt and Legg presented non-uniform UE 
layouts (partly clustered around picocells) according to the 
3 GPP model 4a (T8) . However, this model is not designed 
to adjust the traffic statistically, i.e., the traffic statistics are 
not measured to be used as an input to the traffic generation 
function. 

Dhillon et al., in (19) , proposed a non-uniform UE dis¬ 
tribution model. They start with a higher density of BSs. 
Then they consider a typical UE located at the origin. After 
selecting the serving BS, they condition on this active link and 
independently thin the rest of the BS point process so that the 
resulting density matches the desired density of the actual BSs. 
It should be pointed out that the situations in which UEs are 
clustered, but not necessarily around BSs, are not captured by 
this method. 

The spatial traffic modeling method proposed in (2T), by 
Dongheon Lee et al., suggests that the spatial traffic can be 
approximated by the log-normal or Weibull distribution. This 
paper considers the statistics of the spatial traffic distribution, 
but it does not discuss the modeling of the cross-correlation 
of traffic with BS locations. 

In (23) we proposed new measures for capturing traffic 
characteristics in the space domain. The proposed measures 
can be considered as the analogues of iat in the time domain. 
Thomas point process was used to generate spatial traffic 
patterns with desired characteristics. However, the HetNet 
scenarios are not investigated in (23) . 

In (24) we proposed a novel methodology for the statistical 
modeling of spatial traffic in wireless cellular networks. The 
proposed traffic modeling considered the cross-correlation be¬ 
tween the UEs and BS locations as well as the CoV as defined 
in [23]. The proposed traffic generation method was a density 
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based method with two phases. First, a BS-biased non-uniform 
density function for the entire field was generated, then the 
desired point pattern was produced based on that density 
function. It should be noted, however, that the generation of 
the density function for all the points in the field (as required 
in the method proposed in ||24)) is computationally intensive. 
Moreover, the model proposed in [[24) is not directly applicable 
on HetNets because the density function is calculated for a 
homogeneous macro-only scenario. 

In this paper, we propose a superior traffic generation 
method in comparison to that in |24| which is also applicable 
to HetNet scenarios and study the impact of spatially hetero¬ 
geneous traffic on heterogeneous infrastructure. The proposed 
method is computationally efficient since it does not require 
the generation of a density function. Some of the key results 
in our previous works (23) and (24} are also presented in 
this paper in a comprehensive and coherent way with more 
analytical detail. 


C. Organization of the Paper 

The remainder of this paper is organized as follows: In 
Section [II| our traffic modeling methodology is presented. 
In Section [InJ new traffic measures and appropriate statistics 
are introduced. The traffic generation process is described in 
Section IV The numerical results for traffic modeling and 
network performance analysis are presented in Section |VJ and 


the concluding remarks are made in Section VI 


II. Traffic Modeling Methodology 

Heterogeneous wireless cellular network models receive 
traffic patterns in the time domain as their input for per¬ 
formance analysis. The key performance indicators such as, 
average user data rates, system spectral efficiency, and outage 
probability, are calculated as the system outputs. The traffic 
patterns are usually generated by random processes (e.g., PPP, 
MMPP, HMM, ...), which are called traffic generators (TGs). 

TGs, in turn, receive a number of input parameters (TGIPs) 
to generate traffic patterns with various statistical properties. 
Therefore, a main track of research in traffic modeling has 
been to fit TGIPs to generate traffic patterns with the desired 
statistical properties (e.g., measured statistics from a given real 

traffic trace) ® (25), (26) 

Figure [T] demonstrates our proposed traffic modeling 
methodology in the space domain. A somewhat similar 
methodology has been used in the literature for traffic model¬ 
ing in the time domain GD 

The proposed methodology is described as follows: 

• The desired statistical properties of traffic as the modeling 
inputs are fed to the traffic generator. 

• Since the traffic generator expects TGIPs, the desired 
statistics are first translated to the appropriate TGIPs. 

• Traffic with known characteristics is then generated to be 
used for network performance analysis. 

• The translation look-up table is developed off-line in 
advance: 

- The entire range of the feasible TGIP values are fed 
to the traffic generator and the traffic pattern statistics 


are measured at the output. Therefore, a complete 
map from TGIP values to the corresponding traffic 
statistics is produced. 

- The conversion function from TGIPs to the traffic 
statistics is obtained by fitting methods. 

- The look-up table from traffic statistics to the cor¬ 
responding TGIPs can be prepared by inversion of 
the map or the function from TGIPs to the measured 
statistical properties. 

A mathematical description of the process is given in Section 

IEH 

The method described above is capable of generating a wide 
range of heterogeneity possibilities. However, an appropriate 
next-step is to tune the model according to the real world 
measurements and to extract the accurate model parameters 
for generating realistic traffic patterns. In (27| , we used the 
maps of Paris, France, obtained from OpenStreetMaps (28), 
to study the spatial traffic heterogeneity of outdoor users in 
the denser areas of the city center. 


III. Traffic Measurement 


A perfect model will require all the statistics of the traffic 
measures for the generated traffic patterns and those for 
the real traffic traces, including the cumulative distribution 
functions (CDFs) and auto-correlation functions, to match. 
Therefore, a perfect match is not practical as this requires 
the use of extremely complicated models with a very large 
number of parameters. As a result, simplified models are 
commonly used which only consider the first few moments 
of the traffic. In the time domain, usually the mean, CoV, and 
auto-correlation, and rarely the third moment, are considered 
to match with those of the real traffic trace. 

For the measurement of traffic patterns, first, an appropriate 
measure must be selected to capture the traffic properties. In 
the time domain, iat is the most popular and well accepted 
measure. 

In Section |III-A the equivalent measures for the statistics 
of the spatial traffic are introduced. The spatial traffic statistics 
are described in Section IIII-BI 


A. Traffic Measures 

Packet arrivals in the time domain can be modeled by 
a one-dimensional (ID) point process. A fixed iat between 
packets generates the maximum homogeneity (deterministic 
lattice). An exponentially distributed iat generates complete 
randomness (PPP in ID). For generating sub-Poisson pat¬ 
terns (patterns with more homogeneity than Poisson), one 
way is to generate a perfect lattice, then apply a random 
displacement (perturbation) on its points (29) , (30} . Various 
models for generating super-Poisson patterns (patterns with 
more heterogeneity than Poisson) have been proposed in the 
literature which are mostly based on hierarchical randomness 
and Markov models 0, (n, (nj. 

A ID point pattern in the time domain can be measured 
mathematically in many different ways. One may use the 
interval count, N(a, b\ = Nb—N a (N t representing the number 
of points arrived before time t ), which is a density-based 
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Off-line calculation to generate look-up table: 



Fig. 1. Modeling procedure: The desired statistical properties of traffic as the modeling inputs are translated to the appropriate TGIPs, and traffic with known 
characteristics is generated to be used for network performance analysis (lower dashed box). The look-up table for translation is generated in advance via 
off-line calculations (upper dashed box) [24) . 


measure and divides the whole domain into smaller windows 
and counts the number of pattern points in each window. A 
disadvantage of the density-based measures is that they are 
parameterized by the window size. Finding an appropriate 
window size is itself a challenging question and does not have 
a unique answer for all applications. 

The inter-arrival time, iat, f = T i+ 1 — Ti (T* representing 
the time of arrival of point i), is the most popular and best 
accepted measure because it is distance-based rather than 
density-based, and it considers the distance between every two 
neighboring points in the domain. The CoV for iats is defined 


where p x > 0 and a T are the mean and standard deviation of 
iats, respectively (note that pi > 0 if t > 0). For a perfect ID 
lattice, the constant iat has C T = 0. For a ID Poisson pattern, 
Cj = 1, since for an exponential distribution with parameter 
A, the standard deviation and the mean are both p I = cr I = A. 
Sub-Poisson processes have 0 < C T < 1 and super-Poisson 
processes have C T > 1. 

The UE locations in a heterogeneous wireless cellular net¬ 
work in space domain can be modeled by a two-dimensional 
(2D) or three-dimensional (3D) point process. A very inclusive 
review of Point processes in space domain is conducted in 
m- A fixed distance between points generates perfect homo¬ 
geneity (deterministic lattice). On the contrary, the Poisson 
distribution generates complete randomness. For generating 
sub-Poisson patterns, one way is to generate a perfect lattice, 
then apply a random perturbation on its points g9)-(3T). 

As mentioned above, in time domain, the distance-based 
measure iat captures heterogeneity by one non-parameterized 
real value C T . In multi dimensions, however, there is no natural 
ordering of the points, so finding the analogue of the iat is not 
straightforward. There are many density-based heterogeneity 
measures in the literature such as Ripley’s K-function and pair 
correlation function ED, but they are all parameterized. For 


introducing distance-based measures, there is the problem of 
defining the ’next point’ or the ’neighboring points’ in multi¬ 
dimensional domains. 

The first and simplest candidate for characterizing a neigh¬ 
boring point in a multi-dimensional domain is to consider 
the nearest-neighbor. This leads to the nearest-neighbor dis¬ 
tance measure (32). However, the nearest-neighbor distance 
measure in ID time domain is not the analogue of the iat 
because it considers the min{/^/^ + i} for every point T^. 
It is shown in Section El in our simulation results that the 
nearest-neighbor distance fails to capture the process statistics 
in multi-dimensional domains because this measure ignores 
the neighbors other than the closest one. The next candidate 
is the distance to the k th neighbor. However, determining 
k globally is not possible because every point may have a 
different number of neighbors. 

In the following, novel distance-based and non- 
parameterized measures in the space domain are proposed 
based on the Voronoi and Delaunay tessellations. 

Definition 1. Voronoi Tessellation p3\ p. 1 ]: Given a point 
pattern P = {pi,P 2 , tn d-dimensional space R d , the 

Voronoi tessellation VT = {c Pl , c P2 ,..., c Pn } is the set of cells 
such that every location, y G c Pi , is closer to pi than any other 
point in P . This can he expressed formally as 

c Pi = {y e : \y ~Pi\ < \y-Pj\ for i,j e . (2) 

Definition 2. Delaunay Tessellation p4] p. 11]: The Voronoi 
tessellation in R d has the property that each of its vertices is 
given by the intersection of exactly dpi Voronoi cells. The 
corresponding dpi points define a Delaunay cell. So the two 
tessellations are said to be dual. 

Figure [2] demonstrates a pattern of points with its Voronoi 
tessellation (dashed lines) and Delaunay tessellations (solid 
lines). 

Definition 3. Natural Neighbor / j35] p. 3]: Every two points 
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Fig. 2. Voronoi and Delaunay tessellations: For a point pattern represented 
by bullets, the associated Voronoi tessellations (dashed lines) and Delaunay 
tessellations (solid lines) are illustrated. 


sharing a common edge in Voronoi tessellation or equivalently 
every two connected points in Delaunay tessellation of a point 
process are called ’natural neighbors’. 

Natural neighborship gives an inspiration of neighboring 
relation in multi-dimensional domains and leads us to the ana¬ 
logues of the well accepted iat measure in multi dimensions. 
Various statistical inferences based on different properties of 
cells generated by these tessellations can be considered for the 
measurement of a point pattern. 

The Voronoi cell area or Voronoi cell volume V is the first 
natural choice. For a lattice process, all the cell areas in 2D 
or cell volumes in 3D are equal which results in C v = 0. 
The statistics of the Voronoi cells for a PPP (Poisson-Voronoi 
tessellation) are well investigated in the literature |33|], |36]]- 
f39| . Square-rooted Voronoi cell area in 2D or cube-rooted 
Voronoi cell volume in 3D can also be considered. We are 
interested in Voronoi cell area V because it can be considered 
as an analogue of iat in the time domain. 

The main focus of this paper is the spatial aspect of traffic 
distribution because the temporal aspect is well-investigated 
in the literature. The combination of the temporal domain and 
the spatial domain, which is a very challenging problem, can 
be considered as an extension of this paper. However, it is 
critical to understand the heterogeneity in the space domain 
thoroughly before proceeding to the combined domain. In that 
case, some other metrics such as the inverted Voronoi cell area, 
l/V, could be used to normalize the user target rates in the 
time domain. Therefore, a combined metric such as R/V can 
be of interest in the combined domains; this remains as an 
interesting future work. 

The next proposed measure is the Delaunay edge length 
E. The statistics of Delaunay tessellations is investigated in 
|40|-|[42|. The mean value of the lengths of Delaunay edges 
of every point can also be considered. 

A Delaunay tessellation divides the space into triangles or 
tetrahedrons in 2D and 3D, respectively. The area distribution 
of the triangles or the volume distribution of tetrahedrons can 
determine the properties of the underlying pattern. 

The Voronoi and Delaunay tessellations can be applied on a 
ID process which models traffic in time domain. In this case, 
the introduced distance-based measures are converted to time 
domain measures. The basic statistics of these measures for a 
PPP in one, two, and three dimensions, and their analogues in 
time domain, are summarized in Table [I] 

In order to use the above mentioned measures as an ana¬ 


events events 


▲ 





Fig. 3. Domain similarity: Realizations of processes with sub-Poisson 
(0 < C < 1), Poisson (C = 1), and super-Poisson (C > 1) characteristics, 
respectively, from left to right in time domain (top) and space domain 
(bottom). 


logue of iat, one needs to normalize their CoV with the 
CoV values of iat in the time domain. For the complete 
homogeneity case, the CoV values are already zero, same as 
iat in the time domain. To normalize the CoV values of the 
complete random case to 1, it is required to divide the measure 
by the values presented in Table [I] Figure [3] demonstrates 
realizations of processes with sub-Poisson, Poisson, and super- 
Poisson characteristics. 


B. Traffic Statistics 

Assuming that UEs have the same altitude, spatial traffic is 
modeled as a 2D point pattern U G M 2 which is generated 
by a generator point process T> [/ . In this paper, the mean, 
y rn , and the CoV, C m = a m //i m , are the desired statistics of 
traffic where m is the traffic measure and cr m is the standard 
deviation of m (the third-moment and auto-correlation are 
stated as future extensions in Section [VI|. 

Along with the mean and the CoV, which capture the 
heterogeneity of traffic, a very important statistic of traffic 
in space which affects the network performance is the bias of 
the UE distribution to the BS distribution (i.e., the correlation 
between the two distributions). BS locations can be modeled 
by a superposed 2D point pattern B = Uf=i ^k^k £ 
where K is the number of tiers of BSs in heterogeneous 
infrastructures and B & is the set of BSs of type k operating 
with power P^ and generated by a point process <f> fc with 
density A&. The weighted-Voronoi tessellation of BSs divides 
the entire field into Voronoi cells associated with each BS 
consisting of the area closer, in terms of received signal power, 
to that BS than to any other BS. A sample realization is shown 
in Fig. [4ja). 

To measure the joint distribution of UEs and BSs, we define 
the following potential function: 

Definition 4. Potential Function: Every point (x,y) in the 
field is associated with a potential value P(x,y) G [—1, +1]. 
The P function must have the following properties: 

1) P(x, y) = +1, for cell center points, 

2) P(x,y) = — 1, for Voronoi cell edge points, 

3) If a, p ( x > y) dxd v = °> v *> 
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TABLE I 

Basic statistics of distance-based measures for a PPP in one, two and three dimensions and their analogues in time domain: i is the process point index, A is 
the exponential distribution parameter for inter-arrival time and A is the mean intensity of point processes (23) . 


Distance-based measures 

Time domain 

Statistics 

ID 

2D 

3D 

Nearest-neighbor distance (G) 

min{ Ii,I i+ i} 

Mean (p) 
Variance (a 2 ) 
CoV (C) 

0.5A -i 

0.25A -2 

1 

0.5A _U - 5 

0.0683A -1 

0.653 

0.5539A -U - 33 

0.04A -0 - 66 

0.364 

Voronoi cell area/volume (V) 

A+A + l 

2 

Mean (p) 
Variance (a 2 ) 
CoV (C) 

A" 1 

A" 2 

1 

A” 1 

0.28A -2 

0.529 

A” 1 

0.18A" 2 

0.424 

Delaunay cell edge length (E) 

h 

Mean (p) 
Variance (a 2 ) 
CoV (C) 

A" 1 

A" 2 

1 

1.131A" 0 - 5 

0.31A -1 

0.492 

1.237A -0,33 

0.185A -066 

0.347 



(a) BS weighted-Voronoi tessella- (b) potential distribution 
tion (see Acknowledgment) 


Fig. 4. The weighted-Voronoi tessellation of BSs divides the whole field 
into cells associated with each BS consisting of all points closer, in terms of 
received signal power, to that BS than to any other BS. Every point in the 
field is then associated with a potential value between —1 and +1, according 
to {3}. 



Fig. 5. The d(x , y) is the distance of the point (x, y ) to the associated cell 
center, and D(x,y) is the length of the line connecting the associated cell 
center to the associated Voronoi cell edge through point ( x,y ). 


where i is the index for BSs and Ai is the Voronoi cell area 
associated with BS B P 


The first property assures that if a UE is at the cell center, its 
correlation with BSs is measured as +1. The second property 
assures that if a UE is at the cell edge, its correlation with 
BSs is measured as —1. The third property assures that if all 
UEs are distributed homogeneously (independent from BSs) 
in the area, their mean cross-correlation with BSs is equal to 
0 . 

One may identify many functions which satisfy the require¬ 
ments of the potential function defined above. In this paper, 
we consider the simplest polynomial function. 


Lemma 1 . The simplest polynomial function which satisfies 
the requirements of the potential function is 


P(x,y) 


~2(rf(a:,y)) 2 
(■ D(x,y )) 2 


(3) 


where d(x,y) is the distance of the point (x,y) to the 
associated cell center, and D(x,y) is the length of the line 
connecting the associated cell center to the associated Voronoi 
cell edge through point (x, y) as shown in Fig. [5] 


Proof. See Appendix. 

Figure Qb) illustrates the potential distribution associated 
with the Voronoi tessellation in Fig. [4Ja). 


Using this potential function, the joint moments 

E [P*A J '] , i,j> 0, (4) 

and joint central moments 

E[(P-A* P ) < (A- A t A ) i ],i,j>0, (5) 

can be calculated, where p p is the mean value of P, p A is the 
mean value of UE density A, and E [x] is the expected value 
of x [431. To have a normalized measure, we use correlation 
coefficient which is defined as 


P = 


( 6 ) 


where 


°PA = E [(P - n P )(A - )] 


(7) 


is the covariance of P and A, a p is the standard deviation of 
P, and <t a is the standard deviation of A. So, p can be defined 
as 

= _ - ma )(P(x,y) ~ _ 

JIJIiHx, y) - /u A ) 2 da;<A/)(//(P(:E, y) - fj, P ) 2 dxdy) 

( 8 ) 

For a UE pattern U (a realization of the correlation 
coefficient is calculated as 

Sue*y ^ u 


P = 


\U\ 


(9) 
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where P u is the potential value at point u and \U\ is the 
number of points in U. A pattern with p = +1 means that 
all UEs have gathered at the cell centers, a pattern with p = 0 
means that the UE distribution is independent from the BS 
distribution, and a pattern with p = — 1 means that all UEs 
have gathered at the cell edges. 


IV. Heterogeneous Traffic Generation 


A basic traffic generation method is presented in Section 


IV-A which is further improved in Section IV-B 


A. The Basic Method 

UEs in a heterogeneous wireless cellular network are at¬ 
tracted to social attractors (SA) such as buildings, bus stations, 
shopping centers, and other social places. Furthermore, with 
the proliferation of small-cells and the deployment of BSs 
near potential SAs, the distance between SAs and BSs is de¬ 
creasing. The requirement is a traffic generation method which 
is able to generate most of the possible traffic distributions 
including the following: 

• situations in which UEs are not attracted to SAs at all and 
are completely random, while SAs in turn are distributed 
independently from BSs, 

• situations in which UEs are not attracted to SAs at all 
and are completely random, while SAs are close to BSs, 

• situations in which UEs are highly attracted to SAs, while 
SAs are distributed independently from BSs, 

• situations in which UEs are highly attracted to SAs, and 
SAs are close to BSs. 

Covering the entire range of the above discussed cases, our 
proposed traffic generation method is described as follows. 
Obviously, the proposed method is just one way of generating 
adjustable heterogeneous and BS-correlated traffic; i.e., it is 
not the only way. More accurate methods can be developed 
based on the real world data analysis. 

Starting with a Poissonian distribution of BSs ( B ), including 
macro-BSs, pico-BSs, and femto-BSs, an independent Poisso¬ 
nian distribution of SAs (5), and an independent distribution 
of UEs (U), we move every SA (S 7 ) towards its closest BS, 
in terms of the received signal power ( Bs i ), by a factor of 
a G [0,1]; so the SA’s new location (S^ ew ) is calculated as 

§? ew = aB Si + (l-a)Si. (10) 



Fig. 6. The CDF of the distribution of UEs in the Voronoi cell area of 
SAs: (a) When (3 is deterministic, the UEs are moved towards SAs but the 
UE distribution shape is fixed and the Voronoi cell edge area of SAs remains 
empty with no UEs. (b) When (3 is a random variable with normal distribution, 
even for high (3 values, the Voronoi cell edges are not empty and there is a 
low probability for UEs existing at cell edges. 


B. The Enhanced Method 

An undesired property of the basic method described in 
|IV-A| is that when UEs move towards the SAs, some areas of 
the network become empty having no UEs. In other words, the 
UE distribution shape stays unchanged and it only shrinks to 
a smaller area. Figure [6ja ) shows the CDF of the UE density 
of network after moving UEs towards SAs with different (3 
values using the basic method. Even for small values of (3 
(e.g., /3=0.4), the SA Voronoi cell edge area has no UEs. 

To resolve this undesired feature of the basic method, we 
introduce an enhanced method as follows. Instead of having 
a fixed value for (3, in the enhanced method we model (3 as a 
random variable with 


j8~N(/ jLfrcrp), (12) 

where the mean pp G [0, 1] indicates the closeness factor 
of UEs to SAs (an accurate model based on the real world 
measurements may adopt another distribution for /3, however, 
since we dont have access to real traffic measurements, in this 
paper we choose normal distribution). In ( [l2| ), the a g value 
should have the following characteristics: 

• /3 —y 0 =£’• ap —y 0, 

• (3 —y 1 =r’- <j p —y 0, 

• /5 —^ 0.5 =4> is maximized. 

The <jp value should also be selected in such a way that the 
probability of (3 falling outside [0,1] should be negligible. In 
this paper we set (jp to be 


Then we move every UE (Ui) towards its closest SA, in 
terms of the Euclidean distance (5*y.), by a factor of (3 G [0 1]; 
so the UE’s new location (U™ ew ) is calculated as 

Ur w = f3S Uz + (l-f3)Ui. ( 11 ) 

The initial distribution of UEs can be Poissonian (C = 1) or 
deterministic (C = 0). In this paper, to generate sub-Poisson 
traffic (C < 1), we start with a deterministic lattice. 

The proposed method is general enough to generate traffic 
with negative bias to BSs, i.e., UEs gathering at cell edges. To 
generate this type of traffic, UEs must be moved to cell edges 
instead of cell centers. 


thus the probability of (3 falling outside [0, 1] will be 0.1% 
as shown in Fig. [7] 

Figure [6jb) shows the CDF of the UE density of network 
after moving UEs towards SAs with different /? values using 
the enhanced method. 

Figure [8] illustrates various scenarios with different char¬ 
acteristics which can be generated by the proposed traffic 
generation method. 

In our modeling methodology, the TGIPs are a and (3 which 
are the internal parameters of the traffic generator, and the 
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Fig. 8. Traffic distribution scenarios: By regulating only two parameters, a (the closeness indicator of social attractors to BSs), and (3 (the closeness indicator 
of UEs to social attractors) from 0 to 1, traffic distributions with various properties (CoV and BS-correlation) can be generated. The big circles denote 
macro-BSs, the triangles denote pico-BSs, the lines are the edges of the weighted-Voronoi tessellation of BSs, squares denote social attractors and small 
circles represent UEs. The range of the area is from 0 to 1000 meters in each dimension. 


traffic statistics of interest are C and p. To make the mapping 
(i.e., the look-up table generation) from the statistics (C and 
p in our method) to the TGIPs (a and /3 in our method), we 
first generate traffic patterns with different values of TGIPs 
and then measure the associated statistics. Using a fitting 


procedure, the resulting mapping can be expressed as 


C = F 1 (a,p), 

(14) 

P = F 2 (a,0). 

(15) 
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Fig. 7. (3 is distributed with normal PDF with a mean value of up. 


F\ and F 2 are monotonically non-decreasing functions due 
to the fact that the CoV and correlation coefficient values are 
both non-decreasing with increasing a and [3. Therefore, the 
final maps from the desired statistics to TGIPs can be obtained 
by the corresponding inverse operations as follows: 

a = H 1 (C,p), (16) 


and 

0 = H 2 (C,p). (17) 


V. Numerical Results 

Following the methodology described in Section [II] the 
traffic measurement method discussed in Section [III] and the 
traffic generation model described in Section [IV] the simula¬ 
tion results are presented in this section. The traffic modeling 
results are presented in Section |V-A[ and the heterogeneous 
wireless cellular network performance analysis results are 
presented in Section V-B| 


Fig. 9. To compare the proposed distance-based traffic measures with the 
existing nearest-neighbor distance measure, we fixed the a value to zero (a = 
0) and changed the (3 value from 0 to 1 and calculated the CoV of different 
measures. Nearest-neighbor distance cannot capture the traffic heterogeneity 
and stays constant in super-Poisson region. 



o 

R3 


0.3 0.4 0.5 0.6 0.7 

(3 (UEs movement towards SAs) 


—Delaunay Cell Edges 
-a-Voronoi Cell Areas 
-e-Nearest Neighbor Distance 


Fig. 10. To be consistent with time domain traffic measurement, we divide 
the CoV of the space traffic measures by the convergence value at Poisson 
patterns. The normalized CoV values are used in this paper. 


A. Traffic Modeling Results 

Considering a two-tier superposition of BSs (i.e., macro-BSs 
and pico-BSs), 10 macro-BSs and 20 pico-BSs are distributed 
in a 1000 m x 1000 m square field. The distribution of BSs is 
assumed to be PPR 50 SAs are distributed by an independent 
PPP. UEs are distributed using the traffic generation method 


described in Section IV with parameters a and /3. The simu¬ 
lation is repeated for 1000 random drops. 

The first step in traffic measurement is to select the appro¬ 
priate traffic measure. To compare the proposed spatial traffic 
measures with the existing distance-based nearest-neighbor 
distance measure, with a fixed at 0, we increased (3 from 0 to 
1 and calculated different traffic measures for each f3 value. 
Figure [9] illustrates the CoV values of various distance-based 
traffic measures introduced in this paper. 

The Co Vs of all the measures are normalized to 1 at the 
Poisson case by dividing the CoVs by their expected values 
(presented in Table [I]). Figur e [T0| shows the normalized results. 

As it can be seen in Fig. |10| the nearest-neighbor distance 
measure cannot capture the traffic heterogeneity as it stays 
constant in the super-Poisson region. On the other hand, the 
two proposed measures capture the traffic heterogeneity for 
all CoV values. The Voronoi cell area is preferable to the 


Delaunay cell edge length for two reasons: First, there exists 
a Voronoi cell area associated with each user, while a Delaunay 
cell edge length cannot be associated with one particular UE. 
Secondly, the slope of the Voronoi cell area with respect to (3 
is higher than the slope of the Delaunay cell edge length with 
respect to /?; thus the Voronoi cell area can capture the traffic 
heterogeneity with a higher resolution. In the remainder of this 
paper we use the Voronoi cell area as the traffic measure. 

The next step in traffic modeling is to generate a map from 


TGIPs to traffic statistics. Figures |TT] and [12] demonstrate the 
calculated CoV and correlation coefficient values, respectively, 
for different values of a and /3, both ranging from 0 to 1. 

An important observation from the maps in Figures |TT] and 


12 is that the CoV and BS-correlation are not independent. 


In other words, the range of achievable CoV values for each 
correlation coefficient value is different. This is due to the 
fact that when UEs are biased towards BSs, UEs are attracted 
towards the points of interest, and this automatically shapes 
clusters; which results in a degree of clustering and increases 
the CoV. 

To make a comparison with the existing models in the 
literature, we should note that the prevailing spatial model 
used for UE distribution in the literature is the uniform PPR 
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Fig. 11. The normalized CoV of Voronoi cell area is calculated for traffic 
generated with different values of a and f3. 



Fig. 12. The correlation coefficient (refer to j9j) between UEs and BSs is 
calculated for traffic generated with different values of a and (3. 


To compare the PPP model to our model we should say 
that the PPP corresponds to just one specific point in the 
whole space of the possible spatial distribution situations 
(CoV (heterogeneity) =1 and p (BS-correlation) =0). There 
are few models in the literature which consider heterogeneous 
or BS-correlated traffic modeling. To the best of the authors’ 
knowledge the most appropriate of those is the one presented 
in GD which is a BS-correlated model. However, the model in 
m is capable of generating a limited sub-space of the whole 
traffic possibilities. Figure [T3] illustrates the feasible region of 
CoV and correlation coefficient values generated by our model 
versus the model presented in (T9| . 

The next step in traffic modeling is to translate the desired 
statistics to appropriate TGIPs which can generate traffic 



Fig. 14. The horizontal axis shows the desired normalized CoV. The desired 
normalized CoV is mapped to the associated a and f$ values which are used 
to generate traffic. The vertical axis shows the calculated normalized CoV 
from the generated traffic patterns (for each desired C value, the traffic is 
generated for all feasible p values and the CoV is averaged out). 



Fig. 15. The horizontal axis shows the desired p. The desired p is mapped 
to the associated a and f} values which are used to generate traffic. The 
vertical axis shows the calculated p from the generated traffic patterns (for 
each desired p value, the traffic is generated for all feasible normalized CoV 
values and the correlation coefficient is averaged out). 


with the desired statistics. Towards that end, the statistical 
properties of the generated traffic are measured. Since the 
process of traffic generation is a random process (e.g., /3 is a 
random variable), the measured statistics are not constant and 
have a deviation around their expected values. The deviation 
of the measured statistics from the desired statistics determines 
the accuracy of the method. Figures 14 and 15 illustrate the 
measured statistics versus the desired statistics. 



Fig. 13. The feasible normalized CoV values for different values of 
correlation coefficient are shown. Traffic with high correlation with BSs cannot 
have low normalized CoV values because high correlation means that UEs 
are gathered at cell centers. All the region above the solid line is the possible 
traffic generation using our proposed method. The square-marked line shows 
the possible traffic generation using the model presented in [l9| . 


B. Network Performance Analysis Results 

The performance of a downlink LTE cellular network is 
analyzed in this section. In a field of 1000 m x 1000 m, UEs 
are distributed using the traffic generation method described 
in Sectio n [Tv] The simulation parameters are summarized in 
Table [II] 04J7 (45) . 

The mean value of the rates of all the network users, and 
the coverage probability with a minimum SINR threshold of 
10 dB versus the feasible traffic statistics (for 1000 drops) are 
presented in Fig. [16] and Fig. [IT] respectively. 

The following are the observations made from Figures [l6| 

and [171 

• As the correlation coefficient between UE and BS loca¬ 
tions increases, the mean user rate and coverage proba¬ 
bility also increase as expected, due to the fact that the 
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TABLE II 

Simulation parameters. 


PARAMETER 

VALUE 

Mean number of UEs 

1000 

Number of macro-BSs 

10 

Number of pico-BSs 

20 

Number of SAs 

50 

BSs distribution 

PPP 

BS antenna height 

10 m 

Number of drops 

1000 

Bandwidth (downlink) 

20 MHz 

Noise power per RB 

-174 dBm/Hz 

Carrier frequency 

2.5 GHz 

Total macro-BS transmit 

power 

37 dBm 

Total pico-BS transmit 
power 

17 dBm 

Path loss and shaddowing 

Based on UMi scenario [1] 

BS antenna gain (boresight) 

17 dBi 

UE antenna gain 

0 dBi 

Time domain traffic model 

Full buffer 

Antenna model 

Omni-directional 

BS down tilt 

12 degrees 

UE antenna height 

1.5 m 

Shadowing model 

log-normal with std: LoS: 3, 
NLoS: 6 

Fading model 

no fading 

Scheduling 

Proportional fair [46 ] 



Spatial Traffic correlation coefficient with BSs 


Fig. 16. The mean user rates is shown for the feasible normalized CoV and 
correlation coefficient region. 



Spatial Traffic correlation coefficient with BSs 

Fig. 17. The coverage probability with minimum SINR threshold of 10 dB 
is shown for the feasible normalized CoV and correlation coefficient region. 


received signal strength from the serving BS is increased 
and that from the interfering BSs is decreased. As a result, 
the SINR statistics are improved. 

• For independent UE and BS locations (or, for low correla¬ 
tion coefficient values), as CoV increases, the mean user 
rate and coverage probability decrease. This is due to the 
fact that as CoV increases, the UEs become concentrated 
in some localities, and as a result, some BSs don’t 
have any UEs connected to them, i.e., the resources of 
these BSs are wasted. On the other hand, some BSs are 
overloaded with more than the average number of UEs 
which have to share the available resources. 

• For correlated UE and BS locations (high correlation 
coefficient values), as CoV increases, the mean user rate 
and coverage probability also increase. In these scenarios, 
with increasing CoV, UEs move more and more towards 
the BSs (cell centers); as a result, the SINR statistics are 
improved. 

VI. Concluding Remarks 

In this paper, a new statistical approach for modeling the 
spatial traffic in heterogeneous cellular networks is introduced. 
A number of novel distance-based traffic measures in space are 
proposed which can be considered as the equivalents of the iat 
in the time domain. Only two statistical parameters are used to 
regulate the heterogeneity and bias of traffic towards the BSs. 
Stochastic geometry is used to generate realistic and adjustable 
traffic. The effects of the realistic traffic modeling on the 
performance of heterogeneous wireless cellular networks is 
illustrated. 

This work can be extended in a number of directions. It is 
important to measure the auto-correlation of traffic in space 
as a fundamental statistic which affects the performance. It 
is also beneficial to consider other point processes and other 
Potential functions (refer to ([3])) for traffic generation. Various 
techniques, such as the user-in-the-loop (UIL) (47|-||49), can 
be used to improve network performance for the heterogeneous 
and BS-correlated traffic. Finally, a joint time domain and 
space domain modeling may be developed to obtain a more 
comprehensive traffic model. 

The traffic modeling method presented in this paper is 
general enough to be applied in other wireless networks as 
well. For instance, the relaying and routing problems in ad hoc 
and Wi-Fi networks are closely related, and are dependent on 
the distribution of the network terminals. In wireless networks, 
the cell switch-off framework is also highly dependent on the 
spatial distribution of the traffic. 

It is interesting to note that various characteristics of the 
user distributions in wireless networks can be observed in 
other phenomena as well. For instance, the similar problems of 
heterogeneity of distribution or self-similarity arise in micro- 
cosmic level investigations in chemistry and particle physics 
and in macro-cosmic level studies in astrophysics. 

Acknowledgment 

The authors would like to express their gratitude to Dr. 
Harpreet Dhillon (Virginia Tech, USA) for providing MAT- 
LAB code for generating the weighted-Voronoi tessellation, 



















































12 



Fig. 18. Every BS Voronoi cell can be divided to sub-triangles which 
are each comprised of the BS as a vertex and the lines connecting the 



and L is the length of Voronoi cell edge and 6 is the angle 
associated with the BS vertex. Increasing the number of sub¬ 
triangles, 6 tends to 0 and cos 6 tends to 1. So, 




( 20 ) 


Now, assume that P is a polynomial function. With substitu¬ 
tion of P function in the conditions, we can see that P cannot 
have a degree of one. So, the minimum degree for P is two. 
With substitution of a second-degree polynomial function in 
the potential function conditions, the P function is calculated 
as 

—2x 2 

~^ + 1 ’ (2D 

which can be generalized to any arbitrary point (x,y). This 
proves that the simplest polynomial P function is 


P ( X y ) = zMeM. + 1 

{X ' V) (D(x,y )) 2 + ' 


( 22 ) 
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Fig. 19. The potential value is assumed to be constant on each line 
parallel to the cell edges. So, the P value is only a function of x. 
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Appendix 

Proof of the Simplest Polynomial Potential 
Function 

First, note that every BS Voronoi cell can be divided to 
sub-triangles which are each comprised of the BS as a vertex 
and the lines connecting the BS to the edge. An example is 
shown in Fig. [TS] Increasing the number of sub-triangles to 
large numbers assures that the whole area is covered. 

Also note that for macro-BSs, some parts of the Voronoi 
cell might be covered by picocells. However, since the sum 
of potential value inside each picocell is 0 and the pico-BSs 
are distributed uniformly, it does not affect the sum potential 
values in macrocells. 

If the rules of potential function hold inside each sub¬ 
triangle, we can infer that they hold for entire Voronoi cell 
because the potential value of each cell is simply the sum 
of potential values of sub-triangles. Therefore, the required 
condition for potential function is as follows: 

• P{x,y) = + 1 for cell center points, 

• P(x,y) = — 1 for Voronoi cell edge points, 

• If A, F ( x ’ v) dxd v = 

where A i are the sub-triangles. 

Since the potential value is assumed to be constant on each 
line parallel to the cell edges (as shown in Fig. [T9]), it means 
that the last condition can be written as follows: 

pX PCX pX 

/ / P(x,y)dxdy = / cxP(x)dx , (18) 

Jx=oJy=0 J x =0 

where constant c is 

C=^-7i, (19) 
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